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Abstract: We investigate the polarization properties of high harmonic 
generation from polar molecules with a linearly polarized field. It is found 
that elliptically polarized harmonics are observed in a wide spectral range 
from the plateau to the cutoff. Further analyses show that the nonsymmetric 
structure of the highest occupied molecular orbital is the origin of ellipticity 
of the harmonics. The results provide a method for generation of large- 
ellipticity XUV pulses, which will benefit the application of HHG as a tool 
of detection in materials and biology science. 
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1. Introduction 

High harmonic generation (HHG) has attracted significant attention for producing coherent 
attosecond pulses in the XUV regime |[l}[5] and detecting the molecular structure as well as 
attosecond electronic dynamics ||6}|8| . The power spectrum of high harmonic emission has been 
extensively studied in previous experimental and theoretical investigations ||9}|TT|. Recently, 
the polarization characteristics of the harmonic emission are also actively studied due to the 
potential applications of elliptically polarized harmonics. It is found that information about the 
molecular system participating in the harmonic generation can be imprinted in the polarization 



characteristics of harmonics 1 12 - 14 1. Therefore elliptically polarized harmonics can provide a 
sensitive probe of the molecular system. For example, Y. Mairesse et al. [15 1 both theoretically 
and experimentally demonstrated that elliptically polarized harmonics can be applied in HHG 
spectroscopy of oriented N2 to track multichannel dynamics during strong-field ionization. In 



addition, the search for elUpticity in HHG is also motivated by the generation of elliptically 
(even circularly) polarized attosecond pulses in the xuv regime p6| . 

Several attempts have been made to obtain large ellipticity in the harmonics. It has been 
shown that the nonlinearly polarized harmonics can be obtained using an elliptically polarized 
driving laser |17|. In such case, the harmonic efficiency drops significantly with increasing el- 
lipticity of the driving laser fTE\. The elliptically polarized harmonics are also observed when 
the oriented molecule is driven by a linearly polarized laser | [T9] - |23| . Previous works have 
proposed several origins of the nonzero ellipticity of the harmonics generated from oriented 
molecules. For example, both two-center |'20','2T1 and channel (16) interference, which can 
result in structural and dynamical minima respectively, are responsible for the elliptically po- 
larized harmonics. In detail, the structural and dynamical minima in the parallel component of 
the harmonics make the amplitude of the parallel and perpendicular components comparable. 
Meanwhile, a phase jump by ^ tt of the parallel component is accompanied by the interfer- 
ence minimum position, leading to the phase difference between the two components passing 
through 7z/2. Therefore the harmonics in the vicinity of the interference minima are elliptically 
polarized. Recently, Ramakrishna et al. jT?] demonstrated that the molecular potential also 
plays an important role in determining the polarization properties of the harmonics. In previous 
investigations, the targets are either atoms or non-polar molecules such as H2, CO2, O2, and 
N2. However, the polarization characteristics of the harmonics generated from polar molecules 
are seldom studied to the best of our knowledge. The orbital symmetry is one of the basic prop- 
erties of molecules and can embody the difference between the polar and non-polar molecules. 
Nonetheless, its influence on the elliptical polarization is also seldom investigated. 

In this paper, we investigate the polarization characteristics of HHG from the polar molecules 
with a linearly polarized laser and analyze the influence of the structure of the highest occupied 
molecular orbital (HOMO). Our results show that for oriented polar molecules strongly ellip- 
tically polarized harmonics are observed in a spectral range from the plateau to the cutoff at 
the orientation angles ranging from 30° to 80°. Further analyses reveal that the nonsymmetric 
structure of the HOMO is responsible for those polarization characteristics of the harmonics in 
the case of oriented polar molecules. 



2. Theoretical model 

In this paper, we use the example of CO to investigate the polarization characteristics of the 
harmonic emission from oriented polar molecules driven by a linearly polarized laser Our 
coordinate system, depicted in Fig. 1(a), defines the polarization vector of the driving pulse 
as the space-fixed z axis, with the molecular axis zl lying in the y-z plane. The orientation 
angle between the molecular axis (zl) and the driving laser polarization axis (z) is denoted as 
G. Note that we only need to consider the y, z components of the harmonic emission, since 
the X component along the propagation direction of the driving laser cannot be phase matched 



1 19 22 1. In Fig. 1(b), an illustration of the HHG ellipse is presented. The major, minor axes of 
the ellipse are denoted as a, b respectively. represents the rotation angle of the major axis (a) 
of the HHG ellipse with respect to the polarization axis (z) of the driving laser field. Ellipticity 
of HHG £ is defined as the ratio of the length between the minor and major axes, i.e. e = b/a. 

The simulations are carried out using the extended strong field approximation (SFA) model 
| [24) . With its analytical and fully quantum-mechanical formulations for harmonics, one can 
explicitly and directly analyze the influence of the orbital symmetry of the HOMO. In present 
work, the main focus is the ellipticity of the harmonics, which is dominantly influenced by 
the electronic dynamics. The rotational dynamics are not taken into account by assuming a 
perfect orientation, which is analogous to |20 21]. More sophisticated theory that takes both 



the rotations and the molecular potential into account is presented in recent works y4 23 25 1. 



Fig. 1. (Color online) (a) An illustration of the laboratory coordinate system. (x,y,z) is 
the laboratory frame,the driving laser field propagates along the x axis and oscillates along 
the z axis. The orientation angle 6 is positive for clockwise rotation from the z direction, 
(b) An illustration of HHG ellipse. The a, b axes are the major, minor axes of the ellipse 
respectively. <p is the rotation angle of the a axis with respect to the z axis. 



Within the SFA model, the y, z components of the laser-induced dipole moments responsible 
for harmonic emission are expressed as (in atomic units) 



ry{t;e) 
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In these equations, E{t) refers to the electric field of the driving laser pulse, A(f) is its associated 
vector potential, ^ is a positive constant, t' and f correspond to the moments of electronic 
ionization and recombination respectively, psr and Sst are the stationary momentum and the 
quasi-classical action, which are given by 



p.At',t) = -^ ['A{t")dt", 
t — t Jt' 

sAt',t) = {t-t%-^-pi{t',t){t-t')+y',A\nd, 



(3) 



(4) 



where is the ionization energy of the target molecule. g{t;9) represents the ground state 
amplitude at the time t, and can be expressed as p6l[27|: 



g{t;d)^exp 



w{t';e)dt' 



(5) 



w(f';0) is the ionization rate, which is calculated by Molecular Ammosov-Delone-Krainov 
(MO-ADK) model for aligned molecules |28|. d{p) is the field-free dipole matrix element for 



transition from the ground state to the continuum state, where p stands for the momentum of the 
electron, dy^d^ correspond to its y, z components, respectively. Within the single active electron 
approximation (SAE), this transition dipole moment is given by |13j 

d{p;e)^{^Q{x,yx,0)VWp)- (6) 

Here y]io(x,y,z\ 0) represents the ground state of the target molecule, i.e. the highest occupied 
molecular orbital (HOMO). The HOMO is obtained using the Gaussian 03 ab initio code | [29) 
with the 3-21G basis set. 1//^ = exp {ipz) refers to the electronic continuum state with a momen- 
tum p. The harmonic spectrum is then obtained by Fourier transforming the time-dependent 
dipole acceleration d{t; 9). The equations for the two components are: 



T 



ay{q;d) — — j ay{t;B)exp{—iqO}t)dt, (7) 



1 



T 



a,{q;e) = - a,{t-e)exp{~iq(Ot)dt, (8) 
1 Jo 

where a(f; 9) = r(t; 9), T and CD are the duration and frequency of the driving pulse, q corre- 
sponds to the harmonic order. The spectral intensity ( / ) and phase ((j?) of harmonic components 
are given by Ij{q) = \aj{q)\^ , (Pj{q) = aig[aj{q)] , j — y,z- The ellipticity £ is determined by 
the amplitude ratio and the phase difference of the two components pT] : 



/l+r2-Vl+2r2cos25 + r4 

where r = ^lyjh and 5 = (Py — (p^ . The range of the ellipticity is < £ < 1 . The linear, elliptical, 
and circular polarization correspond to£ = 0,0<£<l, and £ = 1 respectively. If 5 is or 
TT, the HHG will be linearly polarized. For other values of 5 between and n, the HHG will 
be elliptically polarized. In order to attain high ellipticity, the phase difference must be around 
7t/2. 

3. Result and discussion 

We study the ellipticity of HHG from oriented polar CO molecules driven by a linearly polar- 
ized laser.A 30-fs 800-nm driving pulse with an intensity of 2 x 10^^ W /cm^ is used in our 
simulation. The electric field of the laser pulse is expressed as 

E{t)^Eocos{cot)sin^{y% (10) 

where £0 is the amplitude of the driving field. Figure. 2(a)presents the ellipticity of the harmonic 
emission at two orientation angles 30° (red line) and 45° (black line). One can see that for the 
two orientation angles, strongly elliptically polarized harmonics are observed in a wide spectral 
range covering the plateau to the cutoff. In Fig. 2(b), the phase difference (blue line) and the 
ellipticity (black line) of the recombination dipole moment 1 13 30) for CO oriented at 45° 



are also presented. Here the phase difference is shifted into [— TT, 0]. Interestingly, the relative 
phase between the y, z components of the dipole moment is close to —n/2 for all the electronic 
momenta. As a result, the ellipticity of the dipole moment is quite large, as shown by the black 
curve. Furthermore, it is also shown that the ellipticity curve of the dipole moment in panel 
(b) is similar to that of the harmonics at the same orientation angle (black line) in panel (a). 
The electronic momentum is transformed into the harmonic order by the dispersion relation 




Fig. 2. (Color online) (a) The ellipticity of harmonics generated from CO for the orientation 
angle at 30° (the red line) and 45° (the black line). The horizontal axis represents the 
harmonic order, (b) The ellipticity (the black line) and phase difference (the blue line) of 
the recombination dipole moment for the orientation angle at 45°. The phase difference is 
shifted into [— ;r,0]. The electronic momentum is transformed into the harmonic order. 



qhco = Ep +Ip, wherein Ep is the kinetic energy of the free electron with the momentum p. 
This phenomenon will be discussed latter. 

The ellipticity of the harmonics for CO oriented at other angles from 0° to 90° is also stud- 
ied. The map of the harmonic ellipticity is shown in Fig. 3(a) versus the harmonic order and 
the orientation angle. As shown in Fig. 3(a), strongly elliptically polarized harmonics are also 
obtained at other orientation angles ranging from 30° to 80°. In Fig. 3(b), the ellipticity of the 
dipole moment as a function of the orientation angle and the electronic momentum is presented 
as well. Here the electronic momentum is transformed into the harmonic order using the same 
dispersion relation. One can see that the phenomenon observed in Fig. 2 is not specific to the 
case when the molecule is oriented at angle 45°. As shown in Fig. 3, the ellipticity map of the 
harmonics in panel (a) also exhibits a quite similar shape to that of the dipole moment in panel 
(b). To explain this phenomenon, a theoretical analysis is performed based on the three-step 
model 1 24 31 1 for high harmonic generation from oriented molecules. In this model, an elec- 
tron first tunnels from the ground state, then oscillates in the field, finally recombines with the 
parent ion. The recombination step leads to the emission of a high-energy photon with two po- 
larization components, parallel (z) and perpendicular (y) to the laser polarization. Since the first 
two steps are identical for the y, z components of the harmonics, the phase difference and ampli- 
tude ratio between the two components of the harmonics is mainly contributed by those of the 
recombination dipole moment d*{p). Hence the ellipticity curves between harmonics and the 
recombination dipole moment exhibit a similar shape at the same orientation angle, as shown 
in Fig. 2 and Fig. 3. In other words, the recombination dipole moment d*{p) is responsible for 
the large ellipticity of the harmonics over a wide spectral range for CO. 

From the formula d*{p;9) — {exp {ipz)\r\y/o{f", 9)) for the recombination dipole moment, 
one can see that the phase difference arg(t/*) — arg(£/*) and the amplitude ratio \d*/dt\ of the 
recombination dipole moment are directly related to the spatial structure of the HOMO. In the 
following, we analyze the influence of the HOMO structure on the phase difference of d*{p) 
in detail to grasp a deeper insight into the large ellipticity of the harmonics generated from 
the polar CO molecules. For comparison, we also discuss the case of non-polar N2 molecule, 
which has a symmetric HOMO with a similar geometry (Og) to that of CO. Sectional views of 
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Fig. 3. (Color online) (a) Tiie ellipticity of harmonics as a function of the orientation angle 
and the harmonic order. (b)The ellipticity of the recombination dipole moments versus the 
orientation angle and the electronic momentum which is transformed into the harmonic 
order using the dispersion relation q%m = Ep +Ip. Ep is the kinetic energy of the free 
electron with the momentum p, Ip is the ionization energy of the target. For the two panels, 
the target is the polar CO molecule, the orientation angle and the harmonic order correspond 
to the horizontal and vertical axes respectively. 

the HOMO for N2 and CO at the orientation angle 45° are presented in Fig. 4. As shown in 
this figure, the section of the non-polar N2 molecule in panel (a) exhibits a centro-symmetric 
structure. Then the wavefunction satisfies i/Zq (x, y, z; ) = i/o(jc, — y, — z; 0) (for arbitrary value 
of x). However, for the nonsymmetric HOMO of CO this condition is not satisfied, as shown 
in panel (b). With Eq.(6), one can obtain Re[d*] = and Im[d*] ^ for N2, Re\d*^ 7^ and 

Im[d*j] ^ ij ~ y,z) for CO, respectively. As a result the phase difference of d*{p) for N2 
is either or ±7Z, as shown by the blue curve in Fig. 5(b). Whereas for CO, they can take 
other values including ±7z/2 for all the electronic momenta, as demonstrated in Fig. 2(b). The 
phase difference of the harmonics is dominated by that of the recombination dipole moment, 
as illustrated above. Therefore large ellipticity of the harmonics is obtained in a wide spectral 
spectrum, as presented in Fig. 2(a) by the black curve. Finally, it can be concluded from the 
above discussion that the nonsymmetric structure of the HOMO is the major origin of the 
ellipticity observed in the case of the oriented polar CO molecules. 

Figure 5(a) shows the ellipticity of the harmonics generated from N2 at the same orienta- 
tion angles with Fig. 2(a). In contrast to CO, the nonzero ellipticity of the harmonics are only 
obtained in the vicinity of the structure minima in the HHG spectra, where the positions of 
the minima coincide with those in |30|. This nonzero ellipticity of the harmonics is well ex- 
plained by the two-center interference effects, which is analogous to |16 20 21 j for non-polar 
molecules. Comparison between the cases of CO and N2 also shows that the ellipticity of the 
harmonics of CO is larger than that of N2 for most of the orientation angles, as demonstrated in 
Fig. 2 and Fig. 5 for the orientation angles at 30° and 45° . Therefore our results suggest a poten- 
tial application of HHG from polar CO molecules in generating the broadband large-ellipticity 
harmonics, which will benefit the application of HHG as a tool of detection in biology and 
materials science in the future. 
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Fig. 4. (Color online) A view of the section parallel to the yz plane for N2 in the left panel 
and CO in the right panel. The two target molecules are oriented at 45°. 
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Fig. 5. (Color online) (a) The ellipticity of harmonics generated from the oriented N2 at 30° 
(the red line) and 45° (the black line). The horizontal axis is the harmonic order, (b) The 
ellipticity (the black line) and phase difference (the blue line) of the recombination dipole 
moment for N2 oriented at 45°. The phase difference is shifted into [— ;r,0]. The electronic 
momentum is transformed into the harmonic order. 

4. Conclusion 

We use the example of CO to investigate the polarization characteristics of HHG from the ori- 
ented polar molecules driven by a linearly polarized laser. It is found that the harmonics over 
a wide spectral range are strongly elliptically polarized for the orientation angles ranging from 
30° to 80°. The spectral range of large-ellipticity harmonics is covered from the plateau to the 
cutoff, while for non-polar molecules the harmonics only in the vicinity of the interference 
minimum are elliptically polarized. Further analyses show that the nonsymmetric structure of 
the HOMO is the major origin of the ellipticity of the harmonics observed in oriented CO 
molecules. Since the nonsymmetric structure is a general characteristic of the HOMO for po- 
lar molecules, the conclusions made from the results of CO can be extended to other polar 
molecules. Thus our investigations provide a method for generation of large-ellipticity XUV 



pulses, which will benefit the application of HHG as a tool of detection in materials and biol- 
ogy science. 
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